Quantal interferometry with dissipative internal motion 
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In presence of dissipation, quantal states may acquire complex-valued phase eflects. We suggest a 
notion of dissipative interferometry that accommodates this complex- valued structure and that may 
serve as a tool for analyzing the effect of certain kinds of external influences on quantal interference. 
The concept of mixed-state phase and concomitant gauge invariance is extended to dissipative 
internal motion. The resulting complex-valued mixed-state interference effects lead to well-known 
results in the unitary limit and in the case of dissipative motion of pure quantal states. Dissipative 
interferometry is applied to fault-tolerant geometric quantum computation. 
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I. INTRODUCTION 

Garrison and Wright 0] were first to address the effect 
of dissipation on quantal interference, with particular em- 
phasis on the geometric phase in cyclic motion. By mod- 
eling the dissipative quantal motion with non-Hermitian 
Hamiltonians, they arrived at a complex-valued phase 
concept, the geometric interpretation of which was for- 
mulated in a biorthonormal description |^. This result 
has triggered further work on complex-valued geometric 
phase effects S S S • 

In this paper, we revisit, from the point of view of in- 
terferometry, the concept of phase in the presence of dis- 
sipation. Our point of departure is the concept of gauge 
invariance applied to interference in dissipative systems. 
We focus on a proper treatment of gauge invariance of 
the non-Hermitian description of dissipative interferom- 
etry. When the interfering system carries some internal 
degrees of freedom, consideration of this gauge symme- 
try is shown to lead to a complex- valued geometric phase 
effect for arbitrary input states. This effect constitutes 
the non-Hermitian generalization of the mixed-state ge- 
ometric phase put forward in Ref. 9| and, in the partic- 
ular case of pure cyclic internal states, it reduces to the 
Garrison- Wright phase 0]. 

Geometric quantum computation, first proposed in 
Ref. [^ and experimentally demonstrated in Ref. [ll|, 
has attracted considerable interest recently due to its pre- 
dicted resilience to certain kinds of errors. This attrac- 
tive feature has been analyzed from different perspectives 
in the Abelian case, such as random unitary perturba- 
tions ^21, [l3| and decoherence [MJ (for similar analyses of 
non- Abelian geometric quantum computation, see Refs. 
[iSl ll6l lH IS El). We demonstrate fault-tolerance 
with respect to dissipative decay for an Abelian geomet- 
ric phase shift gate. 

The biorthonormal approach to dissiptative quantal 
motion is described in the next section. In particular, 
we put forward an extension to the mixed-state case, in 



order to pave the way for the subsequent analysis in Sec. 
Ill of dissipative interferometry with internal degrees of 
freedom. Section IV contains an analysis of a nonadi- 
abatic one-qubit geometric phase shift gate undergoing 
dissipative decay. The paper ends with the conclusions. 



II. BIORTHONORMAL APPROACH 

In presence of dissipation, we expect the norm of state 
vectors to change in time. To illustrate this, consider a 
dissipative system modeled by a partial absorber, charac- 
terized by the transmission probability < T < 1. When 
passing through the absorber, any |-0) in Hilbert space Ti 
transforms as |V') -^ v2^|V')j which has norm reduced by 
a factor T. On the other hand, the space of pure states 
is P{n) = n/lC - {0}) [mi, C - {0} being the set of 
nonzero complex numbers, i.e., ji/') and vT|V') should be 
regarded as the same state. 

To deal with the projective structure -P(7i), one defines 
pure states as generalized one dimensional projectors Q 
that take the form 



P = |a)(/3|, 



(1) 



such that TrP = (/3|a) == 1. Here |a), \(3) g H are said to 
be binormalized. Dissipation may be modeled by a time- 
dependent non-Hermitian Hamiltonian H{t), which we 
assume to have a nondegenerate complex- valued discrete 
spectrum. In terms of the linear operators L{t) and R(t), 
being solutions of (?i = 1 from now on) 

iL{t) ^ H{t)L{t), 

iR{t) = H\t)R{t) (2) 

with L(0) — -R(O) ~ I, pure states evolve as 

P^P{t)=L{t)PR\t), (3) 

which is trace preserving for arbitrary P provided that 

R\t)L{t)^I. (4) 

Indeed, we obtain from Eq. Q that 

d 
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dt 



(R\t)L{t)) = 0, 



(5) 



which together with L(0) = i?(0) = / imphes Eq. igj. 
In the example with a partial absorber discussed above, 
we may put L = \fT and R = 1/\/T so that |V')(''/'I ~^ 
L\'4') {iIj\R'^ = |'0)('0|, in concurrence with the projective 
structure Piji.)- 

In many physical situations pure states do not provide 
an accurate state description and one has to resort to 
mixed states. To deal with these cases, we introduce the 
generalized density operator 



N 

E 

fc=i 



Wk\ak){l3k\ 



(6) 



where Wk > are real-valued and sum upto unity so 
that Trp = 1 by requiring (/3fc|a/) — Ski UM- The set 
{\ak),\f3k);k = l,...iV} with J2k=i \o^k){(3k\ = / is said 
to be a biorthonormal complete basis of the N dimen- 
sional Hilbert space Ti, such that \ak) and \Pk) are eigen- 
vectors of p and p^ , respectively, p evolves as 



p^p{t)^L{t)pR\t), 



(7) 



which in conjunction with Eq. Q assures preserved 
trace. It further follows that {\ak{t)) = L{t)\ak)} and 
{|/3fc(i)) = R{t)\Pk)} are nonorthogonal sets of eigen- 
vectors of p{t) and p^'it), respectively, both with time- 
independent semi-positive eigenvalues Wk- 
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|0) 



|0> 



FIG. 1: Interferometer setup illustrating the effect of partial 
absorption and phase shift. 1 — T is the absorption probability 
and X a variable U(l) phase. 



Let z = y/Te'-^ and identify L = z and R = 1/z* so 
that |^)(?/'| ^ L|?/')(i/'l-R^ = IV')(V'I for any |?/') € Hs- 
In other words, z does not affect the local motion in the 
|1) beam, but it may show up in interference. To verify 
this latter point, we note that the effect of z may be 
represented by the operators 



L = |0)(0|+z|l)(l|, 

R = |o)(o| + ^|i)(i| 



(9) 



III. DISSIPATIVE INTERFEROMETRY 

Consider a single beam of particles incident on the 
standard Mach-Zehnder interferometer shown in Fig. ^ 
At each equal-time slice in the interferometer, |0),|1) 
span the Hilbert space TCs describing the spatial beam- 
pair. All horizontal (vertical) beams are denoted |0) (|1)). 
Let 1 — T, 0<T<1, be the absorption probability of 
a static partial absorber in the |1) beam added to it a 
variable U(l) phase x- What is the output intensity con- 
sistent with the biorthonormal description? 

To appreciate what follows, we may first recall that the 
standard analysis |23| yields the output intensity 



(8) 



T 1 ^ 2VT 



in the |0) output channel. It is perhaps tempting to inter- 
pret this result in terms of an interference pattern char- 
acterized by visibility v — 2\/T/{l+T) and phase shift x- 
On the other hand, as noted in the previous section, the 
state along the |1) beam inside the interferometer is unaf- 
fected by the application of the U(l) phase, as well as of 
the partial absorber. This gauge symmetry suggests that 
a mere multiplication of ^/Te^^ in the |1) beam should 
only change the phase in the output, while the visibility 
should remain unity. We now show that the biorthonor- 
mal formalism accommodates a natural notion of phase 
and visibility, adapted to this intuition. 



By analyzing the interferometer, we obtain that the input 
state |0)(0| transforms as 

^ Hcr^ZH|0)(0|Ht^tcr^Ht 

= i(2+i + z)|0)(0| 

+i(2-i-z)|l)(l| 
4 z 



-interference terms 



(10) 



with (Ta; and H the standard Pauli-X and Hadamard oper- 
ator, respectively, acting on Ti^. In the |0) output chan- 
nel, we obtain the z-dependent complex-valued intensity 



X[o] oc 2 H \- z 



(11) 



This intensity displays a singularity at the origin in the 
complex z plane. Physically, this singular point at T = 
corresponds to vanishing interference or, equivalently, 
a situation where the path of the particles is perfectly 
known. The complex- valued phase shift cf) and visibility 
V are defined as 



X[o] oc 1 -I- 1/ cos 0, 



which yields 



o^'t' 



(12) 



(13) 



as desired. 

To further analyze the relation between T\n\ and the 
experimental parameters x ^^id T, let us focus on the 
complex- valued interference term cos (j) = j7ro] • By in- 
troducing the polar decomposition j7[o] = l>7[o]|e~*'', we 
obtain 



tan^ 



1 -T 
1 + T 



tanx, 



\Jm\ = xlcos^x' ^^ ^^' 



AT 



(14) 



which shows that jTro] rotates in the complex plane with 
angular frequency that tends to x in the singular T ^ 
limit where iJ/jojl becomes infinite. In the unitary limit 
T — + 1, it follows directly from Eq. (|12|l that j7[o] oscil- 
lates along the real axis according to the expected cos x- 
Next, we extend the above setup and assume that the 
particles carry some internal degrees of freedom prepared 
in an input state that is described by the generalized den- 
sity operator p, so that the incoming state is character- 
ized by the generalized density operator g^^ = |0)(0| (^ p 
acting on the full Hilbert space TCs ® H. Suppose that 

L = \0){0\^L + z\l){l\(g)I, 

R = |0)(0|®i?+— |l)(l|(g)/ (15) 

with L,R being solutions of Eq. 101 for H{t) ^ H\t), 
are applied between the first beam-splitter and the mirror 
pair, see Fig. |21 and the complex- valued z must in general 
contain absorption in order to fully exploit the complex- 
valued structure of the phase shift resulting from the non- 
Hermitian transformation of the internal motion. The 
output state becomes 



gout = UBUAfLUBgiuUTjRTUT,,UTj 



(16) 



with XJm — (Jx® I and Ub = H (g) /, where a^ and H as 
above acting on Tis- The operators Um, U^, L, and R 
act on the full Hilbert space TCs^'H. L and R correspond 
to the application of L and R along the |0) path and the 
z operation similarly along |1). Direct evaluation yields 
the output state 



Pout = ^|0)(0|® (Lpi?t + p+-Lp + zpi?t 

+ i|l)(l|®(Lpi?t 
+interference terms. 



p Lp ~ zpR^ 

z 



In the |0) output channel, we thus obtain 

X[o] oc 2 -h -Tt[Lp] + zTr[pR^ . 



(17) 



(18) 



|1) 




FIG. 2: Interferometer setup illustrating the interference ef- 
fect due to dissipative internal motion. The dissipative mo- 
tion is modeled by the non-Hermitian Hamiltonian H{t) that 
acts on the internal degrees of freedom of the interfering par- 
ticles. 2; is a nonzero complex number. 



with z — e 
yields 



20 



7^ 0, e C. Comparing Eqs. lfTS|l and ifT^ . 



„i$ 



V 



/ Tr [Lp] 



/TT[Lp]Tv[pRt].^ 



(20) 



which is the desired complex- valued generalization of rel- 
ative phase and interference visibility put forward in Rcf. 
[9l|. In the case where H{t) is Hermitian, we have that 
L = R = U is unitary and we obtain 



„z$ 



V = 



lTr[Up] 
Tr[pf/t]' 



/Tt[Up]Tt[pW], 



(21) 



which is consistent with Ref. [PI in the case where p 
is Hermitian. Note, however, that $ may be complex- 
valued in the general unitary case where |/3fc) ^ \ak)- 

We now introduce a concept of geometric phase asso- 
ciated with the path C : i e [0,t] -> p{t) = L{t)pR^{t). 
The basic observation for this purpose is that there is 
an equivalence set A of operator pairs L{t),R{t) that all 
generate C, namely those of the form 

N 



L{t) = L{t)Y,Zk{t)\ak){Pk\ 



k=l 

N 



R{t) = R{t)Y,47TM{au\ 

fe=l ^kK^I 



(22) 



In analogy with Eq. H12|) . we may now introduce the 
complex-valued relative phase $ and visibility V as 



I[o] oc 1 + V cos 



$ 



(19) 



where Zk{t), t e [0, r], are nonvanishing complex numbers 
such that 2:fc(0) = 1. The equivalence set A is the proper 
non-Hermitian generalization of that in the unitary case 
introduced in Ref. [23 • The existence of A expresses the 



gauge symmetry of dissipative motion of mixed quantal 
states. We may identify {Lll(t), i?ll(i)} G A fulfilling 



(/3fc|i?llt(i)Lll(i)|«fc>-0, k^l,...,N. 



(23) 



These constitute parallel transport conditions in a fiber 
bundle with structure group (C — {0})^. Substituting 
L^t),R^t) = L{t),R{t), with L{t),R{t) given by Eq. 
(j^ . into Eq. ^^, we obtain 

Jr^\ ^ g-X"{/3fcl«*(t)i(t)l"fc>rft ^ g~/„"(/3fc(t)|afc(*)>dt (24) 

where we have used that zj| (0) = 1 . Putting this into the 
relative phase, we identify the mixed-state generalization 
of the complex-valued geometric phase factor as 



Tr[Lll(T)p 



Tr[pi?llt(r)] 






(25) 



We may verify that e 



»ri 



'Lit)Mt) = ^' \L{t)Mt) for any 
L{t),R(t) e yl. Thus, e*r is a property of C. 

Let us now analyse two important special cases. First, 
assume that H{t) is Hermitian, so that L{t) = R{t) = 
U{t) is a one-parameter family of unitarities. In this 
case, Eq. H25|l takes the form 



e^r ^ 






(26) 

where \ak{t)) = U{t)\ak) and |/3fc(i)) = U{t)\pk)- Equa- 
tion (|26|l is consistent with the mixed-state geometric 
phase in unitary evolution 9] when p is Hermitian, but 
r may still be complex- valued in the general unitary case 
where |/3fc) ^ \ak). 

Secondly, in the pure cyclic case, defined by p = |q;)(/3|, 



L(r)|a) = e<|a) and i?(T)|/3) = e< ,^ 
complex number [l|, we obtain 

In analogy with Ref. |2J| , we may put 

|5(i)) = e-'/W|a(i)), 
m)) = e-'/-W|/3(i)) 

such that /(r) - /(O) = C- This yields 

(m\a{t))dt, 
which is consistent with Ref. [l|. 



with (^ some 

(27) 

(28) 
(29) 



IV. FAULT-TOLERANT GEOMETRIC 
QUANTUM COMPUTATION 

Abelian nonadiabatic geometric quantum computation 
has been proposed j25i | in order to achieve high-speed 
fault-tolerant implementations of quantum gates. Here, 
we analyze the resilience of a nonadiabatic one-qubit ge- 
ometric phase shift gate to dissipation. 

The physical scenario for the dissipative phase shift 
gate is an unstable two-level atom interacting with an 
external electric field that rotates uniformly around the 
z axis. By using the rotating wave approximation, the 
internal quantal motion of the atom may, in the rotating 
frame, be described by the non-Hermitian Hamiltonian 



H = i^{-n-il)^z = 2^^z, 



(30) 



where rj is the detuning, 7 is the average decay rate, and 
for simplicity we have neglected small off-diagonal terms 
in the weak-amplitude limit of the transverse electric field 
in the x — y plane. Furthermore, CTz = \g){g\ — |e)(e|, 
where \g) and |e) are the unperturbed ground and first 
excited atomic state, respectively. Solving Eq. Q yields 



L{t) = e 
R{t) 



—iujazt/2 



.t/2 



(31) 



Assume 



p=l±L\a+){a+\ + ^\a^){a^\, (32) 

where < r < 1 is the nonzero length of the Bloch vector 
of the qubit and 



|a+) = cos-|.g)+sin-|e). 



l«-> 



-sin-|g) +cos-|e> 



(33) 



make an angle 9 ^ H with respect to the z axis. In 
absence of decay, i.e., 7 = 0, a geometric phase shift gate 
of the form 



U\ 



J{l~ag)2Tr{l~cos9) 



(34) 



where ae — |a+)(a+| — |q;_)(q;_|, is obtained after one 
period r = 27r/?7, by eliminating the dynamical phase, 
e.g., using refocusing technique. The dependence upon 
the solid angle 27r(l — cos 61) expresses the geometric na- 
ture of C/[6l]. 

Now, we may notice that \a±) undergoes noncyclic evo- 
lution for alH > in the presence of decay, except in the 
trivial case where 9 — 0. Yet, we may still compute the 
noncyclic complex- valued phases by using Eqs. (|20|l and 

First, let ip = lot be the complex- valued total preces- 
sion angle. By using Eq. (|20|l . we may compute the 



relative phase and visibility as 
$ = — arctan 
V = ./cos^ 



r cos 6 tan — 
2 



f 



2 cos^ sin 



2 (^ 



(35) 



which are in general complex-valued unless the average 
decay rate 7 vanishes. 

Next, by using Eq. H25|l . we may compute the geomet- 
ric phase associated with the path t e [0, r] — > L{t)pR\t) 
as 



r = — arctan 



r tan — 
2 



where 



J7 = 2 arctan 



cos y tan 



V' 



— (^ cos 6 



(36) 



(37) 



is a complex- valued analog of the geodesically closed solid 
angle on the Bloch sphere, appearing in the Hermitian 
case. For small 7T, we may Taylor expand U, and obtain 



Q, = 2 arctan 



cos 9 tan 



TjT 



— TjT COS 9 



-IJT COS I 



sin^ ^ sin^ 9 
(l - sin2 HI sin2 ( 



0[(7r)2](38) 



whose real part is exactly the desired solid angle to second 
order in "/t. This expression further entails that the one- 
qubit gate U[9] in Eq. H34|) is fault-tolerant in 7, in that 
f2 and, since r is independent of 77 and 7, thereby also 
r are robust to second order in the decay for r = 27T/ri, 



I.e., 



n 



27r(l - cos 6*) -I- O 



27r7y 



(39) 



This feature further supports the predicted resilience of 
geometric quantum computation to unwanted external 
influences. A related result in the context of adiabatic 



evolution of open quantum systems has been found in 
Ref. [11 . 

Furthermore, for small jt we may also expand the rel- 
ative phase as 



$ 



arctan 

."fT 



r cos 9 tan 



r]T 



r cost 



2 1 — (1 — r^ cos^ 9) sm 

21 



2 TIT_ 

2 



+0[(7r) 



(40) 



Thus, the imaginary part of $ vanishes for 9 = -k 12 but is 
nonzero otherwise for all 7T > 0. In other words, the rela- 
tive phase, which contains both dynamical and geometric 
contributions, is not fault-tolerant to the dissipative de- 
cay. This suggests that the above demonstrated resilience 
to the present form of dissipative decay is a consequence 
of the geometric nature of F, and that a quantum gate 
based upon $ would be sensitive to this particular kind 
of error. 

V. CONCLUSIONS 

Two-beam interferometry with particles carrying in- 
ternal degrees of freedom has been analyzed by using 
a biorthonormal dynamical description. It leads to the 
notion of dissipative interferometry that may serve as a 
tool for analyzing the effect of certain kinds of external 
influences on quantal interference. In dissipative interfer- 
ometry, phases and visibilities become complex- valued. 
Gauge invariance, adapted to internal mixed input states, 
is defined in terms of path invariance in state space, and 
has been shown to lead to a natural concept of geometric 
phase that generalizes Ref. _9] to dissipative motion as 
well as Ref. _1] to the mixed-state case. Fault-tolerance 
of Abelian nonadiabatic geometric quantum computation 
is demonstrated for a one-qubit phase shift gate, mod- 
eled by an unstable two-level atomic system. It would be 
pertinent to extend the present framework to the non- 
Abelian case in order to deal with the robustness of uni- 
versal sets of quantum gates implemented by geometric 
means. 
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